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possibility that some bosonic components of the field have a finite vacuum ex- 
pectation value. We derive an exact hierarchy of flow equations for the one-line 
irreducible vertices and the vacuum expectation value of the field. We apply our 
method to an interacting Fermi system where the interaction can be decoupled 
in the zero-sound channel and is then mediated by a collective bosonic field. The 
vacuum expectation value of the zero-frequency and zero-momentum component 
of the bosonic field is then closely related to the fermionic density. This can be 
exploited to calculate the compressibility of the interacting system. By using a 
cutoff in the bosonic propagator, the RG can be set up such that the self-consistent 
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rections to this approximation are generated as the remaining degrees of freedom 
are successively eliminated by the RG procedure. 
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1. Introduction 

Recently, we have proposed (in collaboration with L. Bartosch) a new formulation of 
the functional renormalization group approach to interacting Fermi systems which is 
based on the explicit introduction of collective bosonic fields representing the relevant 
collective fluctuations via a suitable Hubbard-Stratonovich transformation 1 . In 
contrast to the conventional renormalization group (RG) approach to quantum critical 
phenomena pioneered by Hertz , in our approach the Fermi fields are not completely 
integrated out, but the RG flow of the coupled Fermi-Bose theory is derived using 
the functional RG equation for the generating functional of the one-line irreducible 
vertices 0]. A similar strategy has recently been adopted by Baier, Bick and 
Wetterich 5 . However, on the technical level our formulation differs from that of Ref. 
[5]. In particular, we succeeded in introducing the RG cutoff parameter A in such a 
way that the asymptotic Ward identities underlying the exact solubility of the one- 
dimensional Tomonaga-Luttinger model are preserved. This has enabled us to solve 
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the infinite hierarchy of flow equations of the coupled Fermi-Bose theory and reproduce 
the exact solution for the single-particle Green function of the Tomonaga-Luttinger 
model entirely within the framework of the RG. 

As pointed out in Ref. PP, our method should be especially convenient to study 
spontaneous symmetry breaking in interacting Fermi systems provided that the order 
parameter can be defined in terms of the vacuum expectation value of a suitable 
bosonic field. However, the exact RG flow equations in the form given in Ref. pQ 
are only valid as long as the fields do not have finite vacuum expectation values. 
We present here a generalization of the RG flow equations derived in Ref. 1 which 
explicitly allows for finite vacuum expectation values of one or several components 
of the fields. Note that this does not necessarily mean that some symmetry is 
spontaneously broken (see Sec.0J. 

Previously, several authors have studied spontaneous symmetry breaking in 
interacting Fermi systems within the framework of the functional RG [51 E| |S]. 
However, in these works a purely fermionic version of the functional RG has been 
used. In view of the bosonic nature of the order parameter field in interacting Fermi 
systems, a formulation where the flow of the vacuum expectation value of the order 
parameter field appears explicitly in the exact hierarchy of flow equations seems to 
be more advantageous, in particular in systems where the fluctuations of the order 
parameter field are controlled by a non-Gaussian fixed-point. 



2. Exact RG flow equations with finite vacuum expectation values 

In this section we consider a general Euclidean action S[&] = 5q[$] + S\ [<&] involving 
some multi-component field $ Q which can have fermionic and bosonic components. 
The index a is a super-label for all quantities that are necessary to specify the fields, 
such as frequency, momentum, spin, and field-type. We write the quadratic part So[$] 
of the action as 



(1) 



where Go is a matrix in all indices and the symbol J denotes integration over the 
continuous components and summation over the discrete components of the super- 
index a. At this point it is not necessary to specify the interaction part Si[$], 
which can contain terms describing interactions between more than two particles and 
couplings between the fermionic and bosonic components of the super-field $ Q . The 
generating functional Q C \J] of the connected Green functions can be defined as the 
following functional integral, 

e QAJ] = J_ f D § e -5 -S 1 + (J,$) ( (2) 

Zo J 

where J a are super-sources associated with the super-fields, and we use the compact 
notation 

(J,$) = / J Q $ Q . (3) 

J a 

The partition function Z$ of the non- interacting system can be written as the Gaussian 
integral 

Z Q = Jd<$> e - So . (4) 
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Suppose now that at least one of the components of the super-field $ Q has a finite 
expectation value (vacuum expectation value) even for vanishing source fields J, 

sg c [j] 



SJ n 



= (*„>|^o=*a^0 



(5) 



We now derive an exact hierarchy for RG flow equations for the one-line irreducible 
vertices and for the vacuum expectation value. The exact hierarchy obtained in is 
only applicable if none of the field components has a finite vacuum expectation value. 
To generalize the hierarchy of flow equations, we proceed as usual and introduce a 
cutoff A into the theory by modifying Go such that the long-wavelength and low- 
energy modes are suppressed. At this point it is not necessary to specify the precise 
cutoff procedure. Differentiating Eq. (J2J with respect to A we obtain pQ 

T\ 



p. n 1 (5Q C i.SGc 




5J6J 



8a In Z 



where we have defined the following matrix in the super-indices, 



SJSJ 



(6) 



(7) 



SJ a SJ, 

To discuss symmetry breaking, it is more convenient to consider the Legendre 
transform of Q c [J] , 

£[4] = (J[*],*)-&[J[S]] (8) 
which has an extremum at $ a = <I>° , 



= . (9) 

Actually, the field $ in Eqs. (JSJ and denotes the expectation value of the field <& 
in Eq. J3J. To simplify the notation, we use the same symbol for both quantities; it 
is understood that from now on we redefine ($<*) — * $«. In order to derive the exact 
RG flow equations of the one-line irreducible vertices r^"*... )C(n in the presence of a 
finite vacuum expectation value $ , we consider the flow equation for the generating 
functional 



r[$] = £[*] + -(A$, Go X A$) 



(10) 



where A$ Q 



6T[$] 



Note that by construction is also extremal at $ a = $^ 



5<t> c 







The one-line irreducible vertices T. 



(n) 



(11) 

which depend implicitly on the vacuum 
expectation value <&°, are then defined via the functional Taylor expansion around < f )0 
(see, e.g., [!|| and references therein), 



Using Eqs. JE 
9 A r = -~Tr 



r[*] = £ 

n=0 

P|l and SGc/SJ, 



5J5J 
(A^^G^ ]) 



L ai,...,a n a i "•' ^^"n ' 

$ Q , we obtain the exact flow equation 
+ d A In Z 



(12) 



1 



(13) 
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Introducing the matrices X, U, G and G via 

S = 



U = 





T 








5<f>5<f> 












T 




T 




<5$<5$ 




<5<M$ 


4>=<[>o 





G-^G^-S, 

G = — G9a[Gq ^G = [1 - GqS]- 1 (9 a G ) [1 - SCop 1 
we may rewrite Eq. (|13(l as follows, 



(14) 

(15) 

(16) 
(17) 



5* 



1. 



-Tr 



ZG T U T {1 - G T U T pj - (A^^aIGo 1 * ]) 
ZGjlF {1 - G^y 1 ] i(4> , [9 A G *]*°) 



(18) 



where [Z] Q| g = <5a,/3Ca with £ Q = ±1 depending on whether <I> Q is a bosonic or fermionic 
field. To obtain the flow equations for the vertices, we should take into account that 
$° depends on A, so that the derivative of Eq. (|12fl with respect to A yields 

9^ = J2-J ■ ■ ■ I d A T£_ an A<S> ai ..... A$ Q „ 

n=0 Jq i Ja " 

-E^t/ /•••/ ri"i 1 ).., Q „(a A $°)A$ Q1 .....A$ Q „. (19) 

n= Q n - J a Jot! J a n 

Equating the right-hand sides of Eq. l|18[l and 1|19[) we obtain the flow equations for the 
one-line irreducible vertices. The vertex r' ' with no external legs is the interaction 
correction to the grand canonical potential. Its flow equation is 



9 A r 



(0) 



-Tr 



zg^ {i - g^}- 1 ] - i(*o, [d A G^m - (d A $°)r« (20) 



The flow equation for the vertex Fa with one external leg is 



(21) 



and the flow equations for the vertices with n > 2 can be written in closed form 
as follows, 



0ai£>. „, OB = JJd^rtX'U - 2 E E <Wh 



-f ...+mj 



/=1 mi,...,m;-l 

X<S Ql> ...,O ml ;a ml + l ami + m 2 :...i«m l + ...+m,_ l + l <*n { 



Tr 



7 pT r (mi+2) T pTp(m 2 +2) T 

1 ai,...,a mi u A Q mi+ i,...,a mi+r , 

-i(n+2) 



pT-p(m,+2)T 

a roi + ... +m ,_ 1+ i,...,o„ 



where the matrices T ai ,...,a n are defined as follows 

rWri+2) i ; _ p(n+2) 

The operator S symmetrizes the expression in curly brackets with respect to indices 
on different correlation functions, i.e., it generates all permutations of the indices with 



}. (22) 
(23) 
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appropriate signs, counting expressions only once that are generated by permutations 
of indices on the same vertex. More precisely the action of <S is given by (to = ^"V, m.j) 

1 



Sai,...,a mi ;...;a m - mi + i,...,a m {A aii ,,^ am } — „ I Z—t S S n c(^ > ) ^«p(l) , ...,ap(m| ' (^) 

lli TO i- p 

where P denotes a permutation of {1, ... , to} and sgn^ is the sign created by permuting 
field variables according to the permutation P, i.e., 

® am = sgn c (P) $ Qp(1) • . . . • $ Qp(m) . (25) 

We now adjust the flowing vacuum expectation value <I> such that for any value of the 
cutoff A the vertex r„ vanishes. In this case $° can be identified with the flowing 
extremum of the effective potential. The above equations (|2U[1 and (|21[1 can be further 
simplified if we choose the inverse free propagator G^ 1 such that 

G- 1 ^ = , [d A Go V = . (26) 

We assume that this can always be achieved, if necessary by including appropriate 
counter-terms in the Gaussian part of the action. Then the last two terms on the 
right-hand sides of Eas. (|20|l and (|21|) vanish and the latter equation becomes 

in aia d A <S>° a = ^Tr [if'G] , (27) 

which relates the RG flow of the vacuum expectation value <I>° to the irreducible self- 
energy S and the three-legged vertex I^ 3 ). Diagrammatic representations of the flow 
equations for Y^ n ' with n = 1, 2, 3, 4 are shown in Fig. ^ 

3. Flow equations in the symmetry broken phase of the Ising model 

Before applying the general equations derived above to quantum systems involving 
both fermions and bosons, it is instructive to examine these equations in the context 
of classical y> -theory in D dimensions, which describes the critical properties of the D- 
dimensional Ising model. In this case our super-field $ has only a single component ip 
and the super-index a labels momenta k (or, alternatively, points r in D-dimcnsional 
space). Our initial action is S[f] = So[f \ + Si[tp], where 

S<M = \j Go\k)tp- k ip h , (28) 
and the interaction part is 

4 

SiM = ^ / <P-k<Ph + ^ / • ■ • / (27r)^(]T h)ip kl ■■■ip ki . (29) 

Here J fe = yJ2k ~ * I (SL)^ ' ^ or ^ ~^ 00 > wnere V is the volume of the system. 
Furthermore, using a sharp cutoff in momentum space for simplicity, the cutoff- 
dependent free propagator is given by 

9(A<|fc|<A ) 

Go(fe) = — 2 , (30) 

CqK 

where 8(A < |fc| < A ) is unity if A < |fe| < A and vanishes otherwise. For the 
Ising model on a D-dimensional hypercubic lattice with lattice spacing a the bare 
parameters are c = (2D)- 1 and r Q = a~ 2 (T - T C MF )/T C MF , where T C MF is the mean- 
field result for the critical temperature. In the symmetry broken phase the k = 
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Figure 1. Diagrammatic representations of the exact RG flow equations for the 
first few one-line irreducible vertices: (a) rW, (b) r( 2 >, (c) I^ 3 ', (d) rW. Shaded 
circles represent one-line irreducible vertex functions, whereas circles labeled by 
G and G represent the full Green function and the single scale propagator defined 
in Eqs. 1161 and 1171 respectively. The small dotted circle with a cross denotes the 
change d/± "3?° of the vacuum expectation value. The action of the symmetrization 
operator S is defined in Eq. 1241 . The arrow on the circles indicates the order of 
the indices on the corresponding correlation function, which is important to keep 
track of the signs, when Fermi fields are present. 
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Fourier component of the field has a finite expectation value. It is then natural to 
shift 



<p% = (2Tr) D 6(k)M . 



(31) 



For the Ising model, the parameter M has the physical meaning of the magnetization 
per unit volume. Note that GQ 1 (k)(p t k — 0, so that the condition l(2"fi|l is satisfied. 
This is the reason why we have included the quadratic part of the action proportional 
to ro into the definition of the interaction term in Eq. i|29|) . Our initial action then 
takes the form S[<p] = So[Aip} + Si[ip°; Aip] with 



r JL M * + ^M 4 
2 4! 



S 1 [<p ;A<p]=S 1 [<p^cp + A V ,]=V 

+ r o 1) A< A) + i f r o 2) (fc)A^_ fe A^ 

1 Jk 

3 

+ ^/ / [ (^) D S(J2k l )T i 3) (k 1 ,k 2 ,k 3 )A^ kl A^ k2 A^ k3 

Jk t Jk 2 Jk 3 f=1 

+ 4/ •••/ (2^) D ^fc,)r[ ) 4) (fc 1 ,fc 2 ,fc3,fc4)A^ fcl A^ fc2 A^ fc 3A^ fcl , 

^- J ki J fc 4 i=1 



(32) 



where the bare vertices now depend on M, 



r 



r M+^M 3 
6 



ro 



T ( Q 2 \k) 

4 s) (fei,fc2,fc 3 ) =«oM 
r[, 4) (fci, k 2 , k 3 , k^ = u Q . 



—M 2 



(33) 

(34) 

(35) 
(36) 



We now require that Tg 1 ^ = 0, so that the initial value Mq of the order parameter in our 
functional RG is the minimum of the effective potential in the Landau approximation, 



Mn = 



for r > 

y/-§r /u for r < 

Note that for ro < we have 

-.(2)/,_\ u 



wo 



-m 2 , 



(37) 



(38) 



so that in this case the interaction part of our initial action can be written in real 
space as 



SM = 



(39) 



With a sharp cutoff in momentum space we obtain from Eq. (|20|) for the interaction 
correction to the grand canonical potential, 



V 
~2 



/ S(\k\ - A) In 

Jk 



coA 2 + rf (fe) 



c A 2 



The order parameter flow equation (|27|) reduces to 



(5 A M A )rf(0) 



G A (fc)r| 3) (fc,-M). 



(40) 



(41) 
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For a sharp cutoff in momentum space, we have 
and hence 

GA{k) = -^A» + rW(Jb) ' (43) 

From Eq. (|22|) . we obtain for the two-point vertex 

d A r^(k) = l -J^G A (k')Tf{k'-k\k : -k) 

- [ G A {k')G A (k' + k)T^ ) {k,~k~k',k')T { l ] (-k',k + k',-k) 

Jk' 

+(d A M A )rf (fe,-fe,0). (44) 

The above system of flow equations is exact but not closed. In order to obtain a closed 
system for the two-point vertex and the order parameter, it is necessary to express the 
vertices r( 3 ) and r( 4 ) in terms of and M. We now propose a simple polynomial 
truncation which is motivated by the generalized gradient expansion for the effective 
potential ;9. . Guided by our initial conditions 1(3511 , 1(36(1 and 1(38(1 , we truncate the 
above equations by substituting on the right-hand sides 

rf(fc) (45) 

r^ 3) (fc 1 ,fc 2 ,fc 3 ) « u A Al A , (46) 

r^\k u k 2 ,k 3 ,k 4 )^u A . (47) 

If we ignore the momentum-dependencies of the vertices also on the left-hand sides of 
the flow equations, this truncation amounts to approximating 

r A [d= Jd D r^(r)-Mlf, (48) 

which is nothing but the quartic approximation to the effective potential in the 
symmetry broken phase to zeroth order in the gradient approximation In this 
approximation we obtain from Eq. 1(41(1 for the flow of the order parameter, 

d A M\ = -3 / G A (k) , (49) 

Jk 

while Eq. I|44|) reduces to 

d A u A = -3u 2 A [ G A (k)G A (k) . (50) 

Jk 

We shall not further analyze these equations, because they are completely equivalent 
to the RG flow equations discussed by Berges et al. [Sj using the truncated effective 
potential 1)48(1 to zeroth order in the gradient expansion. Whether or not the 
polynomial truncation is justified is a different issue; according to Refs. lllj . 
this truncation fails to describe some important features of the symmetry breaking 
phenomenon. With the help of Eq. (|44|1 one can easily go beyond this approximation. 
For example, we may solve this equation iteratively by substituting the truncation 
in Eqs. I(45|l - ((47|l only on the right-hand sides, retaining the momentum-dependence 
generated by the Green's functions ^1 Ej- From the momentum-dependent self- 
energy we can then obtain an estimate for the anomalous dimension. Having gained 
some confidence in our RG equations, we now focus on a more interesting problem 
involving both bosonic and fermionic degrees of freedom. 
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4. RG approach to interacting Fermions via partial bosonization in the 
zero-sound channel 

In the approach developed in ,1„ tad-pole diagrams of the Hartree-type have not 
been considered due to the assumed overall charge neutrality of the system. For more 
general interactions, these diagrams do contribute and can lead to a non-vanishing 
expectation value of the field describing density fluctuations. In this section, this 
effect is treated with the formalism developed above for symmetry breaking. 

We consider a system with a density-density interaction. After a Hubbard- 
Stratonovich transformation in the zero-sound channel the interaction is mediated 
by a real bosonic field ip and the resulting action reads pQ 

S$,i/>M = So$M + So[<p]+SS, l >l>M > (si) 
with the free parts 

S o 0,V] = -y) / $K*Go*(K)il> Kt , (52) 
SoM = \ {MKrX^R^K*. , (53) 

aa' " 

and the coupling between Fermi and Bose fields 

Sl[W,^]=iV/ / ^K+K.a^KafRa + 777777 [/o^PW^W ■ (54) 
a JkJr 2/jV a,*> 

We use composite frequency-momentum indices K = (ioj,k) and K = (iui^k) which 
include fcrmionic and bosonic momenta and Matsubara frequencies. The integration 
measure is f K = k ancl reduces to J |^ {2n) D ' m * ne zero-temperature and 

thermodynamic limit 0, V —> oo. The propagators are given by 

G ,a(K) = ; , (55) 

FoMX) = Ifk^ 1 - s k,o/(PV))]-^9(A b < \k\ < A fli0 ) , (56) 

with an arbitrary fermionic dispersion ^ka = £fcer — p measured relative to the chemical 
potential p 7 the Fourier transformed interaction , and some average Fermi velocity 
vq which is introduced to give Ap units of momentum. Here, S* = s Af=Ab=0 is the 
self energy of the fully interacting system, which has to be introduced as a counter- 
term in order to scale toward the true Fermi surface |14|. We have allowed for the 
possibility of independent cutoffs on bosonic and fermionic degrees of freedom. 

The Legendre transformation is performed for fermionic as well as bosonic fields, 
such that the resulting functional T[ip,ip,<p] generates one- line irreducible vertices 
that are irreducible in the particle propagator as well as interaction lines |T]. At 
finite density, the bosonic field <p acquires an expectation value and should be shifted 
according to fita ~ <Pr + ^Vko- with 

fRa = - i5 Rfl ( P<y , < P"=Y1 fo a 'p"' ■ (57) 

a 1 

By integrating out the bosonic fields, it is straightforward to see that for Ap = Kg = 0, 
the quantity p & is the density of the interacting electron system in the spin channel 
cr, i.e., p a = $ k {ck a Ckcr) , where Cka is the usual second quantized fermion annihilation 
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operator. The zero frequency and momentum parts of the interaction have not been 
included in F to fulfill the property (|26[1 of the free propagators. The functional 
Taylor expansion 



m^] = ±±j^ [ / .../ /_ .../ 

tTn^TTn \ n -> m - JK a' JK'cr'JKta! J K n cr„ J K 1 S 1 Jb 



n—0 m— 



K' n a' n JK 1 a 1 J K n a n J K±S± J K m S„ 



x $K' 1 + ...+K; i J< 1 + ...+K n +K 1 + ...+K m 

xT (2n,n) {K y^ ^ Ky n ;K 1( 7 U . . . , K n a n ; Kxax, K m a m ) 

X^K[a[ ■ ■■■■ ^K^'^Kiat ■ ■■ ^K n a n ^K 1 a 1 Al fiK m a m > ( 58 ) 

defines the vertex functions Y^ 2n ' m \ 
4-1- Fermi surface cutoff scheme 

In the standard RG approach to interacting Fermi systems |TS], the fermionic degrees 
of freedom are integrated out in successive shells around the (interacting) Fermi 
surface. In our formulation this is achieved by setting As — and by taking only 
A = Ap as a running cutoff. In this case our general flow Eq. (|27|l reduces to the 
following flow equation for the expectation values Cp a : 



E 



JK 



(59) 



where £ = — 1. The fermionic self-energy satisfies the flow equation 

+ [ G a (K')F aa ,(K-K')T^(K'a;Ka;K' - K,a) 

JK' 

x Y {2 ^\Ka;K'a\K - K',a) 

-CV/ G al {K')Y^\Ka,K'a'-K'a\KG), (60) 
(j 1 Jk ' 

where for a sharp cutoff, the propagators are given by 

G a (K) = : , , , (61) 

W - 6ccr - S CT (A) 

r (K , Q(A< \^ + K(^^f)\/vq < A ) 

G a (K) = : -. , , (62) 

F aa ,(K) = [Fo(K)- 1 +n{K)]-^, . (63) 

On the right-hand side of the last equality it is understood, that the O-function in 
Fq should be omitted. The fermionic self-energy and the polarization are defined by 
H a {K) = r( 2 <°)(AV; Ka), and IL aff >(K) = T^ 2 \-Ka, Ka') respectively. 

To explore these flow equations in a simple situation, we will now neglect the flow 
of F^ 2 ' 1 ), r^ 4 ' ) as well as II and set these vertices equal to their initial values. This 
approximation then reduces to an RG version of the Hartree-Fock approximation. No 
frequency dependencies of the self-energy are generated, i.e., T, a (K) — Y. a (k), and we 
can thus define a flowing quasi-particle dispersion = £/- a + £ CT (fe). Using Eqs. (|BTjl 
and Ij59(l . we find that the flow of the density is given by 



dAP* = - I 6(\£ ka + K(iO,kF)\/vo-A)n F (Zt)> 

Jk 



(64) 
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while Eq. I|6U|) determines the flow of the quasi-particle dispersion, 
dACL^Y,fo y (d APa/ )+ [ 6(\£ k ,. + K(iO,kF)\/v -A)fZ. k ,n F (& a ), (65) 

a> Jk ' 

where n F (e) = \<^ c + 1] _1 is the Fermi function, and we have explicitly set f = — 1. 
To further simplify these equations, we use a Hubbard-type interaction with f? a = 
foSa.-a' ■ In this case no extra momentum dependencies of the quasi-particle dispersion 
are generated, so that we may set S a (fe) = —Ap a . Our counter-term is then simply 
E* (iO, k F ) = — liniA^o Afi a , and the flowing quasi-particle dispersion can be written 
as £j£ = £k<? ~ Ap a . If the bare dispersion is actually independent of the spin 
projection a, then also Ap a = Ap and p a — p are spin-independent. In this case Ap 
can be interpreted as the flowing interaction correction to the chemical potential. The 
flow equation l|65[) for the energy dispersion then reduces to 

d A Ap = -f (d AP ) . (66) 

Integrating this equation, we find that, for a given value of the cutoff A, the relation 
between the flowing interaction correction Ap to the chemical potential and the flowing 
density is simply 

An = -f p. (67) 

Substituting this result into Eq. Ij64(l . we obtain a closed flow equation for the RG 
flow of the density p of electrons (per spin projection), 

d\P = -Vq[vq{h - f p* + v A-)n F (f Q (p - p*) + v A) 

+ v (jjl - fop* - voA)n F (fo(p - p*) - v A)} , (68) 

where p* = p\=o is the true density (per spin projection) and 

v {e) = I S{e ka - e) (69) 

Jk 

is the density of states of the non-interacting system. At zero temperature, Eq. (|68|l 
can easily be integrated to yield an implicit relation between the true p* and the 
chemical potential p, 

P* = / deu (e), (70) 

J A 

where A — p — fop* — foAo is the lower band edge. From this expression, we obtain 
the usual Fermi-liquid result for the compressibility |16| . 

dp* vo 



dp 1 + favo ' 

where vo = vo{p — foP*) is the density of states at the true Fermi surface. 



(71) 



4-2. Interaction cutoff flow 

Instead of the standard cutoff procedure for interacting fermions ^E] discussed in the 
previous section, we have proposed in Ref. £Q to use a running cutoff in the bosonic 
sector of the theory. Such a cutoff procedure was also used in the original work by 
Hertz 2 on quantum critical phenomena in interacting electron systems. However, 
in contrast to our approach, Hertz completely eliminated the fcrmionic degrees of 
freedom. In the notation used here, we implement the interaction-cutoff procedure by 
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Figure 2. Initial condition for the pure boson vertices in the momentum transfer 
cutoff scheme. The sum is taken over the m\ permutations of the labels of the 
external legs. 



setting Ap = and by using A = as a running cutoff. The initial condition for the 
flow in this interaction cutoff scheme corresponds to a situation in which interaction 
lines are turned off while particle propagator lines are fully functional. In addition to 
the bare three-leg interaction vertex, the only one-line irreducible diagrams that can be 
drawn are closed loops of fermionic propagators. These loops have to be symmetrized 
with respect to the exchange of external bosonic legs in order to obtain the initial 
condition for the purely bosonic vertices as shown in Fig. [2J A more formal derivation 
of the initial condition is given in the Appendix. The result is given by 

r Ao [^, ij,,ip] = + ^7 XXT'Vo^O*' - Tr{hi[i - iG^]} , (72) 



where the matrix tp of bosonic fields is defined in analogy to Eq. i|A.4|) and 

[G q ]k*,k><,> = 8 K M,5 a ^GoAK) , G ,„(K) = [iw - ^j" 1 • (73) 

The first two terms in Eq. 1)72(1 are nothing but the bare interaction Si defined in 
Eq. 1(54(1 . whereas the last term generates the closed fermion loops. 

When expanded in powers of the field ip, Eq. 1(72(1 contains linear terms. We thus 
need to perform the shift p^a — <P%; + ^PRa with YKo- gi ven i n Eq. 1(57(1 , yielding 

r Ao [*/>, ip, <p] = -Tr{ln[i - iG v Aip}} + ^[/n -1 ]™ Ap 0a Ap 0a , -i^AipoaPa 

aa' a 

+ pM, + A^) - TrIlnlGoG- 1 ]} - ^ £ ff P.P.' , (74) 

aa' 

where we have defined the shifted bare Green function 

G- 1 - Go" 1 - ip° , (75) 
with the matrix elements 

[G^Ka.K'a' = $ K .K> 0~a,a'Go,a{K) , G , a {K) = [iw - Ifco-]" 1 , (76) 

and £ka = £,ka + Yla' fo^ P°' • From Eq. ((74(1 the initial condition for the one-line 
irreducible vertices can be obtained by expanding the right-hand side in powers of 
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A(p and by symmetrizing the resulting expressions with respect to the interchange of 
bosonic fields. We obtain 



r 



(0,0) 



-TtiHGoG- 1 ]} 



2 



-pv{ 



E 



In 



1 



E-^T P°P> 
aa' 



.(0,1) 
A ,cr 

.(2,1) 



Pa 



1 + e-PZ* 

G , a {K) , 



- 2 j: 



fo"' P<rP«'} 



T^(K + K,a;Ka;Ka)=i, 
Tf a Q \Ka;Ka) = X$?(K) = V /<T A 



(77) 

(78) 

(79) 
(80) 



and for m > 2, 



T^iK.a, . . .,K m a) = L^\R X , ...,K m ). 



(82) 



All other one- line irreducible vertices vanish initially. Here, we have defined the 
symmetrized fermion loops as 



■E 



K 



GoAK)GoAK + Kp(i) 



The summation is over all permutations P of {1, . 



,(0,1) 



-... + £p(n-l)). (83) 

, n}. If we now demand that 



0, we obtain the relation 



JK Jk , 



(84) 



which is nothing but a Hartree self-consistency condition. Hence, the initial condition 
in the interaction-cutoff scheme is simply the self-consistent Hartree approximation 
for the density. Note that Eq. (|84|l can also have ferromagnetic solutions with ^ p[. 

By solving the flow equations for the purely bosonic vertices in some appropriate 
truncation, our approach then allows for a systematic inclusion of bosonic fluctuations 
around the self-consistent Hartree approximation. Since the fermionic degrees of 
freedom are not completely integrated out in our approach, it is in principle possible to 
use this solution for the bosonic vertices in the flow of vertices with external fermionic 
legs, e.g. to calculate approximations to the single-particle Green's function. 



5. Summary and outlook 

We have presented a unified way of treating vacuum expectation values of bosonic fields 
in the functional renormalization group. In a simple truncation, our flow equations for 
the Ising model in the symmetry broken phase reduce to well known results obtained 
in the context of the local potential approximation. We apply our formalism to the 
coupled Bose-Fermi theory obtained by decoupling the electron-electron interaction in 
the zero-sound channel. The vacuum expectation value of the zero-frequency and zero- 
momentum mode of the bosonic field is in this case closely related to the fermionic 
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density. The dependence of this vacuum expectation value at the end of the RG 
flow on the initial chemical potential allows us to calculate the compressibility of the 
system. Using a fermionic cutoff and a simple Hartree-type truncation of the flow 
equations, we obtain the known Fermi-liquid result for the compressibility. In the 
interaction-cutoff scheme, a cutoff is introduced solely in the bosonic propagator and 
a fermionic determinant has to be evaluated to determine the initial condition of the 
flow. This initial condition is then given by the self-consistent Hartree approximation 
and the flow equations can be used to obtain corrections to this mean-field theory. 

It is straightforward to generalize the partial bosonization approach to the case 
where some other collective bosonic field develops a finite expectation value. For 
example, with the help of a Hubbard-Stratonovich transformation in the particle- 
particle channel, we may study fluctuation corrections to the BCS approximation for 
a superconductor. Note that within the interaction-cutoff scheme, the self-consistent 
BCS approximation should emerge as the initial condition for the coupled Fermi-Bose 
theory, in analogy to Eq. (|84|l . This should be compared with the purely fermionic 
functional RG approach to symmetry breaking developed by Salmhofer and coauthors 
where the BCS self-consistency condition is only generated in the limit A — > by 
the RG flow starting from an infinitesimal symmetry breaking perturbation. 

The truncation scheme used in SecOlto study symmetry breaking in the classical 
Ising model can be applied to any bosonic sector of quantum systems involving both 
bosons and fermions. Such a procedure, which is equivalent to a truncated local 
potential approximation, might be a convenient starting point to study the feedback 
between bosonic fluctuations on the fermionic single-particle excitations. 

Finally, it should be mentioned that an alternative functional RG approach, which 
also seems to be convenient to study spontaneous symmetry breaking, has recently 
been developed by Dupuis [T2|. His approach is based on the functional RG flow of 
the generating functional of the two-particle irreducible vertex functions. Note that in 
fermionic language the one-line irreducible vertex functions introduced in our partial 
bosonization approach are only partially two-particle irreducible in the sense that 
two-particle irreducible diagrams singled out by the Hubbard-Stratonovich field are 
eliminated in favour of a collective bosonic field. 

We appreciate useful discussions with N. Hasselmann, L. Bartosch and S. 
Ledowski. This work was partially supported by the DFG via the Forschergruppe 
FOR 412. 

Appendix: Initial condition in the interaction-cutoff scheme 

A more formal route to the initial condition l|72|) for T in the interaction-cutoff scheme 
uses an additional generating functional for partially amputated-connected Green's 
functions for which external interaction lines are amputated [18j . 

e Spac[j,.j,J] = — f £)[^ ) ^, ) ^] e -s [->?,V']-SoH-S , i[^>^+J]+(j,V)+('/ i .i) , (A.l) 

Zq J 

By elementary manipulation of the functional integral, we can show 

e a P ac[J,j,J] _ e h(z7< p °Ti) — [ D[$ j ^] e -So$,1>]-S 1 $tf,J]+(jrt)M4>j) . (A.2) 
•Z^i J 

Here, [-Fb]_£ CT R'a' = &R+K' qFo,<t<t'{K') is the matrix of the bare interaction, and 
Z^ = J D[ip, ij^e -30 ^'^ is the partition function for non-interacting particles. For 
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the initial condition of the flow, we have / — > 0. This limit can readily be taken in 
Eq. (|A.2fl . The remaining functional integral is Gaussian and yields 

effij, J] = Tr{ln[i - iG J}} - (J, [1 - iGo^Goj) , (A.3) 

where we use a matrix J of bosonic sources containing the matrix elements 

[J]k<t,K'o' = <W' Jk-K' , a ■ (A. 4) 

The first term in Eq. (|A. 3|) generates the closed loops of fermion propagators in Fig. [21 
when expanded in powers of the bosonic sources. The second term generates diagrams 
that contain a continuous fermionic path linking two external fermionic legs. An 
arbitrary number of external bosonic legs are then directly attached to this line. The 
latter diagrams are not one-line irreducible and will cancel in the expression for T. 
However, before we can perform the Legendre transformation, we first need a relation 
between Q pac and the generating functional Q c for connected Green's functions. This 
is achieved by a shift <p — > cp — J in the integration variables in Eq. (|A.1|I and yields 

G c [J, j, J] = So [J] + 5 P ac [J, j, J] , (A.5) 

where we have defined J — F J. The classical field ip is then given by 

<p = TT = J + Fa—f~-\j,3,J\ = J- (A.6) 
oJ SJ 

In the limit of a vanishing interaction the classical field ip thus becomes identical to 
the source field J. Since Gp^° is quadratic in the fermionic sources, the inversion 
necessary to obtain the sources J and j as a function of the classical fields involves just 
a matrix inversion. The remaining Legendre transformation can then be explicitly 
performed and we obtain Eq. 17211 . 
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